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We show that the quasi-EucHdean sections of various rotating black holes in different 
dimensions possess at least one non-conformal negative mode when thermodynamic insta- 
bilities are expected. The boundary conditions of fixed induced metric correspond to the 
partition function of the grand-canonical ensemble. Indeed, in the asymptotically flat cases, 
we find that a negative mode persists even if the specific heat at constant angular momenta 
is positive, since the stability in this ensemble also requires the positivity of the isothermal 
moment of inertia. We focus in particular on Kerr black holes, on Myers-Perry black holes 
in five and six dimensions, and on the Emparan-Reall black ring solution. We go on further 
to consider the richer case of the asymptotically AdS Kerr black hole in four dimensions, 
where thermodynamic stability is expected for a large enough cosmological constant. The 
results are consistent with previous findings in the non-rotation limit and support the use of 
quasi-Euclidean instantons to construct gravitational partition functions. 



I. INTRODUCTION 

Gravitation is a purely attractive force, which has led to a number of conundrums revolving 
around questions of stability. In classical physics, it seems that the question is settled. Gravitational 
collapse cannot, under a wide range of circumstances, be prevented. Although collapse to form 
a singularity happens, it is believed that these singularities will be isolated from observation by 
horizons. Thus the end-point of gravitational collapse is believed to always result in black holes. 
Classically, many black hole solutions are stable, in particular the ones of astrophysical relevance 

[UEIEIIIEIEIIT!. 

Quantum mechanics changes this. In 1974, Hawking discovered that black holes have a tem- 
perature, the Hawking temperature Th, given by 

K 
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where k is the surface gravity of the black hole j8] (we use natural units, so that G = c = h = kB = l 
throughout). For non-rotating black holes, this gives rise to a new instability. Since 

1 

4M' 

where M is the black hole mass, isolated black holes will radiate and lose energy. This will cause 
them to heat up. Since conservation of energy leads to 

M ~ — 

M2 

we see this is a runaway process. When the black hole reaches zero mass, it is presumed to disappear 
completely. The specific heat of the black hole is negative 

c ' 



87rM2' 

a typical sign of instability [9] . 

It is a sign that the canonical ensemble breaks down for such objects leading to doubts as to 
whether a conventional thermodynamic interpretation is possible. However, if instead one looks at 
the microcanonical ensemble, one discovers that it is well-defined. 

If one includes the possibility that the black holes are rotating with angular momentum J or 
have an electric charge Q, one finds that if 

+ QJ'^M'^ + AQ'^M^ - 3M^ > 

then the specific heat at constant J and Q turns out to be positive. 

One expects that these difficulties will be reflected in the path-integral treatment of gravitation. 
Suppose one tries to calculate the canonical, or grand-canonical, partition function. Then one 
needs to integrate over all physical fields subject to certain boundary conditions. Generally for 



the grand-canonical ensemble, one integrates over quasi-Euclidean configurations (see section II B 
for a more complete description of what this means). So, the field configurations must be periodic 
in imaginary time, with periodicity equal to the inverse temperature, and quasi-periodic in the 
complexified azimuthal angle generated by any conserved angular momenta, or quasi-periodic under 
complexified gauge transformations associated to any conserved charge. 

The gravitational path integral based on the Einstein action is not well-defined because of 
its lack of renormalizability. However, at the semi-classical level, it makes sense as an effective 
field theory, perhaps derived from some more fundamental theory such as string theory. This path 
integral has, at first sight, a big difficulty with stability, as the kinetic energy operator for conformal 
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transformations has the wrong sign. However, it turns out that fluctuations in the path integral of 
such a type are gauge artifacts. Much more serious and interesting is the possibility that the gauge 
invariant parts of the fluctuations contribute with the 'wrong' sign to the partition function. This 
is a sign of instability. For the four-dimensional non-rotating black holes, these negative modes 
have been known for some time |10l [TT] . 

In this paper, we extend our knowledge of this type of instability. In Section llTl we describe the 



formalism required to identify the gauge invariant negative modes. In Section III we describe a 
way of finding gauge invariant deformations of a specified field configuration. These deformations 
are not themselves the negative mode, but since they decrease the Euclidean action, they prove 



that negative modes exist. In Section IV we apply our technique to the four-dimensional Kerr 
solution, five and six-dimensional Myers-Perry metrics, the singly-spinning five-dimensional black 
ring and the four-dimensional Kerr-AdS solution. In every case, we find negative modes, except 
for large Kerr-AdS black holes. In Section |V] we look at the thermodynamics of the black holes or 
rings, and see how it matches up with the existence of our negative modes. We use this to make 
speculations about when the thermodynamic approximation is, or is not, valid. 

II. THE GRAVITATIONAL PATH INTEGRAL 
A. The decomposition theorem 

The path integral of Euclidean quantum gravity. 



Z = J Db]e-^[5], (1) 



is constructed from the action 



^[^] = -T^/ d^^V5(i?-2A)-^ / d''-^xJ^)K-Io. (2) 

The first term is the usual Einstein-Hilbert action and the second is the York-Gibbons-Hawking 
boundary term \\.2\ [13], where K is the trace of the extrinsic curvature on dM. This term is 
required for non-compact manifolds A^, as the ones we will study, so that the boundary condition 
on dM is a fixed induced metric, and not fixed derivatives of the metric normal to dM. 

The term Iq can depend only on ^(^^ ^\ the induced metric on dM, and not on the bulk 
metric gab, so that it can be absorbed into the measure of the path integral. However, since we are 
interested in the partition functions of black holes, it is convenient to choose it so that / = for the 
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background spacetime that the black hole solution approaches asymptotically. For asymptotically 
flat black holes [13], the Einstein-Hilbert term is zero and the action becomes 

d'~'xJ^){K-Ko), (3) 
JdM ^ 

where Kq is the trace of the extrinsic curvature of the flat spacetime matching the black hole metric 
on the boundary dAi at infinity. This subtraction renders the action of the black hole finite. For 
asymptotically AdS black holes |14| [T5] , the boundary terms cancel when the background subtrac- 
tion is performed, but the bulk volume integral diverges and requires an analogous subtraction that 
sets the action of AdS space to zero. An alternative view is that Iq should be seen as a counterterm, 
corresponding to the counterterm of a dual conformal field theory (see |16| [T7\ [TBI [T9l [20] ) . 

The gravitational path integral ([T]) is non-renormalisable but we expect meaningful results in 
an effective field theory approach. A different issue is that the action Q can be made arbitrarily 
negative so that the path integral appears to be always divergent even at tree-level. These problems 
can be addressed in the semiclassical approximation, where the path integral is dealt with by saddle- 
point methods. We consider a saddle-point gab, i-G. a non-singular solution of the equations of 
motion, 

- 2A _ 

Rab = ^ _ r^ ab, (4) 

usually referred to as a gravitational instanton. We then treat as a quantum field hab the small 
perturbations about the saddle-point, 

gab = dab + hab- (5) 

This leads to a perturbative expansion of the action, 

I[g\=I[g\+h[h;g]+0[h'). (6) 

The first order action Ii vanishes since Qab obeys the equations of motion, while the second order 
action I2, which gives the one-loop correction, is the action for the quantum field hab on the 
background geometry gab- 

The effective field theory is valid if the background geometry gab has a curvature nowhere near 
the Planck scale. We can also address the issue of the arbitrarily negative action geometries in 
the path integral, called 'conformal factor problem' since it is the conformal direction in the space 
of metrics that is responsible for the divergence. Perturbatively, this corresponds to trace-like 
perturbations hab which lead to a negative I2- The prescription of [21^ is that the integration 
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contour for those perturbations is imaginary. They can then be seen to be irrelevant and don't 
represent physical instabilities. Of physical interest are the instabilities studied firstly in [10], the 
analysis of which we wish to extend to rotating black holes. 

We follow here the procedure in [22], straightforwardly extended to higher dimensions. We will 
decompose the second order action, applying a standard gauge fixing procedure, and show that 
the unphysical divergent modes do not contribute to the one-loop partition function. 

The partition function is 

^i-ioop = e-^[^] I B[h]{G.F.)e-''^^-'a\ (7) 

where (G.F.) denotes all contributions induced by fixing the gauge in the path integral. Hereafter, 
Qab is relabelled as gab and all metric operations are performed with it. The second order action is 
given by 

" It „, 1 



(8) 



where • denotes the metric contraction of tensors. We have defined 



Kb = Kb - ^9abKc (9) 

and 

{Gh)ab = -V'^VcKb - 2Ra'ffKd, (10) 

where the operator G is related to the Lichnerowicz Laplacian Al hy G = Al — 4A/((i — 2). We 
also define the operations on tensors T 

{ST)b...c = -V'Tab^.c, (11a) 

iaTU...c = y(aTb...c)- (lib) 
The second order action l2[h;g] is invariant for the diffeomorphism transformations 

Kb Kb + VaVb + VbVa = ih + 2aV)ab- (12) 
Following the Feynman-DeWitt-Faddeev-Popov gauge fixing method. 



(G.F.) = {detC)6{Ca[h]-Wa). 



(13) 
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We consider the linear class of gauges 



Cb[h] = ( Kb - \gabh\ ] , (14) 



where (3 is an arbitrary constant, so that the Fadeev-Popov determinant (det C) is given by the 
spectrum of the operator 

{CV)a = -V'VbVa - RabV + (| " l) ^a^bV. (15) 

To study the spectrum, let us consider the Hodge-de Rham decomposition of the gauge vector V 
into harmonic (H), exact (E) and coexact (C) parts, 

V = Vn + VE + Vc. (16) 

This induces a decomposition of the action of C, which we denote by Ch for harmonic vectors, Ce 
for exact vectors and Cc for coexact vectors. 

The harmonic part satisfies dVa = and JVh = 0. We can check that 

^''« = -d^'*' 

The spectrum is positive for A < and zero for A = 0, with multiplicity given by the number 
of linearly independent harmonic vector fields. For A > 0, the background solution satisfying Q 
does not allow for harmonic vector fields if assumed to be compact and orientable [23]. Thus, the 
spectrum of Ch is never negative. 

The exact part is such that Ve = dx, where x is a scalar. We can show that 



spec Ce = spec ( 2 



1 \ 2A 

1 □ 



(18) 



J J d-2_ 

where the operator on the RHS acts on scalars, and □ is the Laplacian. For A < 0, the operator 
is positive for /? > 1, being positive semi-definite for A = 0. For A > 0, the Lichnerowicz-Obata 
theorem tells us that the spectrum of the Laplacian on a compact and orientable manifold satisfying 
Q is bounded from above by —2dK/{{d — l){d — 2)), the saturation of the bound corresponding 
to the sphere [23]. This implies that, for A > 0, the spectrum of Ce is positive for (3 > d. 
The coexact part is such that 8Vc = 0. Hence 

CVc = 25aVc (19) 

and the spectrum of Cc can be shown to be positive semi-definite, 

j A'^x^ [Vc ■ CVc] = 2 j d'^x^ [aVc ■ aVc] > 0, (20) 
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with equality for coexact Killing vectors. 

The Faddeev-Popov determinant contribution to the partition function is then 

detC ~ (detCE)(detC'c), 



(21) 



the tilde denoting that the zero modes have been projected out. The harmonic contribution is not 
explicitly considered because, if it exists (A < 0), it is a positive factor dependent only on A and 
on the dimension of the space of harmonic vector fields, as mentioned above; it will not be relevant 
to our discussion. The contribution from the exact part is fundamental since it will cancel the 
divergent modes of the field hah- 



In order to make the results independent of the arbitrary vector w in the gauge fixing (13 1, 
the 't Hooft method of averaging over gauges is adopted. The arbitrariness is then expressed in 
terms of a constant 7 introduced by the weighting factor of the averaging. The final result will be 
independent of 7, as required. The unconstrained effective action for the perturbations is given by 



If[h-g]=h[h;g] + 



6.''x^C^[h]Ca[h] 



1 

16^ 



7 
327r 

lh-Gh+\{l-^){6hf + l 



1 



Sh ■ dh 



1 



{dhf 



(22) 



where we denote h = h'^^. 

We now decompose the quantum field hab into a traceless-transverse (TT) part, a traceless- 
longitudinal (TL) part, built from a vector r/, and a trace part. 



hab = h^b + Kb + -^dabh, 



1 



with 



Kb = '^ioir])ab + -^gab^V- 



(23) 



(24) 



The constant /3, unspecified in the gauge condition (14), can be chosen so that the trace h and 
the longitudinal vector 77 decouple. This requires 



The effective action becomes 



/3 = 2 1 



d- 27 - 1 



d 



7 



-1 



(25) 



1 



:h' ' ■ Gh' 



arj ■ aAiri -611057]+ 

a 



+ 2(1 — 7) ( 6ari ■ Sar] + '^'^^V " Oi^V ~ ^^'^^ " ) + 



—A lar] . arj - ^-{6v) 



+ -hFh 
2 



(26) 
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where the operator F is given by 



d-2 f d-27-l\^ 1. 



Recalhng the choice of /3 (25), we find that the operator on the RHS of the expression (18) is given 



by 4:dF/{d — 2). The contribution of the ghosts (21 ) can be recast as 



det C ~ (det F) (det Cc). (28) 

For the vector r], as we did for V in the ghost part, we perform a Hodge-de Rham decomposition 
into harmonic, coexact and exact parts, 

11 = ^11 + ^0 + VE, (29) 
respectively. Using ?7e = dx, the result for the effective action is then 



Ih^^ -Gh^^ + lhFh+ 



4 M 
+ ^^3^7Aar/H • ar]n - 7a??c • aCcVc - JZT^,'^^'^ ' ' ^^^^ 



where we defined the operator 



Dab = ^aVb - \9ab^- (31) 



Notice that the Hodge-de Rham decomposition of rj in harmonic, coexact and exact parts gives, 
for /i^ff , a decomposition in 2a??H5 2a7/c and 2Dx, respectively. 

Finally, we can evaluate the Gaussian integrals in the partition function to show the dependence 

Zi_ioop ~ (det C)(det G')~^/2(detF)-i/2(jg^(5^^-i/2(^g^^^-i/2 

~ (detG)-i/2(jgt^^)i/2^ (32) 

Again, the tilde on the operators denotes that the zero modes have been projected out. The 
Gaussian integrals are regularised by (^-function methods [22j. It is understood that the spectrum 
of G here is restricted to traceless-transverse normalisable modes. 

Let us review the treatment of the 'conformal factor problem'. Trace- type perturbations make 
the action ([s]) negative. But a detailed analysis showed that these modes do not contribute to the 
path integral. The two factors (Det-F)~^/^ arising from the Gaussian integrals in h and x cancel 
with the Det-F factor arising from the exact part of the Fadeev-Popov determinant. This makes the 
unphysical character of the divergence obvious, at least in perturbation theory. The conclusion is 
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that the non-positivity of the action ([2]) and the resulting apparent divergence of the gravitational 
path integral are fixed by projecting out this contribution. 

The relevant operators are then Cc and G. For a real metric, the operator Cc is positive 
semi-definite, as we have shown above. Once its zero modes are projected out, it contributes a 
positive factor to the final result. The physical instabilities, identified by imaginary contributions 
to the partition function, are only possible if there are negative eigenvalues of the operator G, 
Gh^^ = Xh^T. This was the problem studied in [10^ for the Schwarzschild black hole. We intend 
to extend this treatment to rotating black holes, which requires addressing the problem of complex 
instantons. 

B. Quasi-Euclidean geometries 

The partition function is usually defined as a Euclidean path integral, a sum over real geome- 
tries for which imaginary time t = it is used. However, while static geometries remain real for 
this analytical continuation, the same does not hold for stationary geometries. In the canonical 
formalism, where 'jij is the metric on a constant time slice, is the lapse function and N^' is the 
shift vector required for rotating spacetimes, we have 

ds^ = N^dT^ + -fij (dx* - iiVMr) {dx^ - iN^dr) . (33) 

These geometries have been called quasi-Euclidean. The question is then whether one should 
analytically continue the shift vector (e.g. through the rotation parameters for a Kerr black hole) 
in order to get a real geometry. This is trivial when one considers the instanton approximation to 
the path integral, because the parameters made imaginary can simply be made real again in the 
final result. But when one goes beyond leading order, as is the case in this paper, and considers 
metrics that do not satisfy the equations of motion but are also included in the sum, the choice 
affects the positivity properties of the second order action and thus the convergence of the path 
integral. 

We share the view of [2l] and |25] that the continuations other than the usual t = it lead to 
unphysical parameters. As those authors point out, the leading order instanton action is real in 
spite of being constructed with a complex metric, and it corresponds to the physical free energy. 
The charges and the horizon locus remain the same as in the Lorentzian case. Studying the 
convergence of the path integral for particular imaginary values of the Kerr rotation parameters, 
for instance, bears no relation to the actual black holes. A further argument can be made based 
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on the black ring case. As opposed to the Kerr geometry, the black ring does not possess a real 
section with imaginary time that is regular, since conical singularities cannot be removed |26| [27] . 
The results in the previous decomposition of the metric were obtained for real Euclidean metrics. 



However, the expression (32) should still hold for an appropriate complex contour of integration. 
This contour is specified in a standard way by the steepest descent method. The relevant eigenvalues 
of Cc and G are now determined with respect to a complex metric, i.e. to the physical Lorentzian 
rotation parameters. The semi-positivity of Cc is no longer obvious and we have nothing more to 
say about it. Still, a negative eigenvalue of G is sufficient to cause problems in the definition of 
the path integral and herald an instability. 

III. THE PROBE PERTURBATION 

As was discussed in the previous section, the Euclidean path integral only depends on the 
spectrum of two operators: Gq and G. The latter acts on traceless-transverse (TT) perturbations 
of the metric and will be the object of our attention in this section. In the Schwarzschild case, 
mostly due to the spherical symmetry of the problem, it was possible to determine the negative 
mode by a straightforward method [10 . However, for solutions such as Kerr, Myers- Perry or the 
black ring, this seems challenging, due to the lack of symmetry of the background geometry. 

The approach that we will adopt here is somehow different. To prove that G possesses at least 
one negative mode we only need to show that a particular TT probe perturbation renders the 
operator negative. In order to visualise this more clearly, pick an arbitrary TT perturbation and 
decompose it in eigenmodes of G, 

</'a6 = 5]an0i?. (34) 

n 

We can now construct the Rayleigh-Ritz functional, given by 

/d^xVff<^-0ce Ep^l ■ ^ ^ 

If a perturbation (pab is found such that I is negative, then it must be the case that at least one 
of the Xn is negative. This reasoning can be used to prove that a given instanton has a negative 
mode, but cannot be used to prove the converse. In fact, if X is positive for a particular (pab^ it 
might be the case that the a„ corresponding to the negative eigenmode is small, or even absent, in 



the expansion (34). 

We also have to check that our particular perturbation lies along the path of steepest descent. 
Since we will consider only perturbations that preserve the symmetries of the background, it suffices 
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to check that the components {G(f))ab are real or imaginary if and only if the components /i"* are, 
i.e. if and only if the components of the metric Qab are. This will indeed be the case for our 
perturbations: the components (r, x*) are imaginary and the others are real. Then, each term of 
the sum (f)"'^{G(j))ab is real and this particular direction in the space of perturbations keeps the phase 
of the integrand functional constant, since the second order action is real. This is the condition for 
the steepest descent path. 

It is the objective of this section to construct a TT perturbation that will render I negative, 
starting with a Killing vector field ka- We will focus on pure Einstein gravity, with A = 0, in an 
arbitrary number of dimensions d, where the background field equations are 

Rab = 0. (36) 

For such a class of spacetimes we can construct a probe Maxwell field satisfying Maxwell's 
equations. Since ka is a Killing vector, it obeys to 

Vafcfe + VbK = 0, (37) 
from which we construct the following two-form components 

Fab = Vakb - Vbka = 2V ah- (38) 

This field strength Fah trivially obeys to the Bianchi identities and also satisfies Maxwell's 
equations, since 

y^pab ^ _2Rb^ka = 0. (39) 

In four dimensions, there exists a TT tensor that can be associated with such a field strength, the 
electromagnetic energy-momentum tensor. 

Tab = F^'Fbc - ^gabPP'^Fpg. (40) 



The energy-momentum tensor defined in Eq. (40) is transverse in any dimension, since that only 



depends on the Bianchi identities associated with Fab and on Eq. (39), but the traceless condition 
is only valid in four dimensions. The strategy is to add a transverse component to Eq. (40), which 
will remove the trace. This can be accomplished by introducing an auxiliary scalar field a, and by 
defining 

(l)ab = Tab- {VaVbCr - QabOa), (41) 
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so that 

V^cl^ab = -Rbay""^ = 0. (42) 

Requiring the tracelessness of (j)ab in Eq. (41) gives 



°^ = = (^)(V^A;'?)(V,A;,), (43) 



For an arbitrary ka, Eq. (43) seems hopeless to invert. However, because ka is a Kihing vector, 



there is a simple particular solution. 



The strategy is now clear: given ka we can construct Fab and a by using Eqs. (38) and (44), 



respectively. These quantities are the only ingredients in the construction of (pab, see Eqs. (40) and 



(41 ). We thus conclude that for each A;" we can associate a TT probe perturbation, in an arbitrary 
number of dimensions. 

We are only interested in perturbations that are normalisable, in the sense that 
f d'^Xy/g(f)"'^(j)ab < +00. The obvious candidates for the Killing vector fields in a black hole back- 
ground are either the time translational Killing vector 9,-, or the azimuthal Killing vector d^, 
both guaranteed in the case of an arbitrary black hole solution in pure d— dimensional Einstein 
gravity |28| [29] . d^f, leads to a non-normalisable perturbation, whereas d-r leads to normalisable 
perturbations. Let us prove the last statement. If the spacetime is asymptotically flat, then 
g^^ ~ 1 + 0{l/r'^~^) and g^^i ~ 0{l/r'^~^), which means that ka — 5a — 0{l/r'^^^) and thus 
F°-b — 0{l/r'^^'^). The last statement implies that y/g<j)'^^(pab — 0{l/r^^'^~'^^), from which we can 
see that this mode is normalisable as long as d > 3. 

We also generalised the construction above to include a cosmological constant. However, we were 
only able to check for the particular cases of d = 4, 5 and 6 that the perturbation was TT, in the 
background of Kerr-AdS and Myers-Perry- AdS in the corresponding dimensions. Unfortunately, 
in d = 5 the perturbation turns out to be non-normalisable, and d = 6 is computationally too 
challenging, so that we will only focus on the Kerr-AdS case. The form of the TT perturbation is 
more involved, 

(^ab = fab - {VaVba - gab^a) - ^^F^a '^b)c, (45) 



where 



rab = Fa'Fbc-^-fFP''Fp,, a=^-(^^yk''ka, (46a) 
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Fab = ^ah - ^bka, Fab = ^ ah " V bK, (46b) 

and 

ka = ka — ka, ka = Qabk^ ■ (46c) 

In the expressions above, ka is a Killing vector of the original metric Qab and Qab is a reference 
metric, obtained from the original metric by setting the mass of the black hole to zero, that is, the 
AdS metric. Note that ka is a Killing vector of gab-, but ka is not a Killing vector field of either 
metrics gab or gab- Moreover, Fab satisfies the sourceless Maxwell's equations. The authors strongly 
believe that the perturbation described above should be TT for all Kerr-Schild spacetimes in any 
dimension, whose reference metric gab is that of a maximally symmetric spacetime. 

In the next Sections, we will be able to identify negative modes using this probe perturbation. 



IV. NEGATIVE MODES OF GRAVITATIONAL INSTANTONS 

We will now apply the method described in the previous section to the Kerr black hole |30] . 
the Myers-Perry black hole in five and six dimensions [31], and the five dimensional black ring of 
Emparan and Reall |32] . All asymptotically flat black holes that we have studied have at least one 
normalisable negative mode, which suggests it may be a universal feature. In the last subsection, 
we study the Kerr- AdS black hole. 



A. Kerr black hole 



The complexified version of the Kerr metric is given by 



ds^ = ^ — + ^ 



S2 



ia{r'^ + 0? — A)dT 

7 



+ ^dr2 + sW, 



(47) 



where 



TP' =r^ + cos^ I 



(48a) 



A 



ror -|- a 



(48b) 



and 



p = [r -|- a j — Aa sm d- 



(48c) 
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The Kerr metric written in this way is already in the canonical (ADM) form. Here, ro is a mass 
scale, and is related to the black hole mass by tq = 2M. Black holes require a < ro/2, where 
the inequality is saturated in the extremal limit. The complementary limit corresponds to naked 



singularities. The avoidance of a conical singularity at = ro/2 + y^Jr^j4)~^^a? requires the 
coordinate identification (r, 0) = (r + /3, — where f3 = (r^ + a^)/[27r(r^ — ro/2)] and 

O = a/{r'^ + a^) are the black hole inverse temperature and angular velocity, respectively. 

For this particular case, due to its simplicity, we will present the explicit expression of the TT 
perturbation that we used to prove that Z is negative, 

,2, 



2S4 



1 + 



2a2 sin^ 6 








2ia{r'^ + a^) sin" 









1 














-1 





2ia 








/ 2a2 sin^ 9 





(49) 



This perturbation is clearly traceless and can be checked to be transverse and normalisable. The 



Rayleigh-Ritz functional defined in Eq. ( 35 ) , evaluated for the perturbation ( 49 ) , is 
50rga2 + 15r^)a - 15r^[2a'^ + (2r_^ - ro)2](2r^ - ro) arctan (f- 



2a?r'^ 



(50) 



(4a2 - 3r^)a + 3r§(2r^ - ro) arctan yf 

This expression is always negative and finite, as can be seen in Fig. [T] We conclude that the 
Kerr instanton is unstable for non-conformal perturbations. It is fortunate that this particular 
probe perturbation was able to identify the negative mode. Note that X(o = 0) = — 5/(7rQ) ~ 
-O.Tlro 2 > -0.76ro 2, the negative eigenmode of the Lichnerowicz operator found in [10] . 

B. Five dimensional Myers-Perry 

We will use the complexified form of the line element of the five dimensional Myers-Perry 
solution in the coordinates introduced in 1331. 



ds" = {x + y) — + 



+ - 



xy 



'-[df-ixyd4)2-i{x-y)d(f)i]'^, 



(51) 



\4X AY J y{x + y) x{x + y) 

where the coordinates (f, 0, ip) are related to the canonically defined coordinates (r, (j), ip) via 

r = f -I- i(af -I- a^^i -I- iafa2'^2, 4>i = ai(j)i + aia24>2 and 02 = 0201 + 020^02- (^2) 
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Kerr 




FIG. 1: For the Kerr instanton, X is negative, decreasing monotonically away from a = and evaluating to 
— 12 r'^'^ at extremality \a\ — rQ/2. 



Also, X and Y are quadratic polynomials whose coefficients depend on the mass 
rotation parameters (01,02), 

X = {x + al){x + al) -rl and Y = -{y - al){y - al). 



2M and 



(53) 



Again, to avoid a conical singularity at the horizon, located at the larger real root of X, 
one must require the coordinate identification {T,(l)i,(p2) = {t + f3,(j)i — if3i^i,(j)2 — 'i'P^2), where 
P = 2TT{x^+ai){x_^^+a2)/[y/x^X'{x_^_)] and ilj = ai/{x^+af) are the black hole inverse temperature 
and angular velocities, respectively. To compute the final integral we further need to specify 
the range of the coordinates on our manifold: < t < (3, < (pi, (p2 < 27r, x > x_^ and 
min(af , 02) < y < max(a^, ag). 

The final result for I is cumbersome and we will only graphically represent, in Fig. Q, the 
variation of X as a function of oi/ro and 02/^0. The similarities with the Kerr black hole are 



evident. The shaded area is the parameter space where Eq. (51 1 represents a black hole, that 
is of + + 2|aia2| < Tq, where equality corresponds to extremality, and it coincides with the 
region where I is negative. In contrast with the Kerr geometry, the five dimensional Myers-Perry 
geometry admits solutions with singular (infinite negative) I, the corners in Fig. These are 
the extremal solutions with rotation in a single plane, for which the horizon area is zero and there 
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J To - Five Dimensional Myers-Perry 
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FIG. 2: For the five dimensional Myers-Perry instanton, X is always negative. The dashed lines represent 
contours of constant Xr\. 

is a naked singularity. The leading behaviour near the singular points is well approximated by 

274 



(54) 



609(ai ±ro)2' 

In certain limits, X considerably simplifies, as for instance in the five dimensional Tangherlini 
solution (ai = a2 = 0) [3l], in which case 

24 



19r2' 



(55) 



or in the singly-spinning Myers-Perry black hole limit (02 = 0), 

8(160ai^ - 1040rgaP _^ 2794r^aJ° - 3927rgaf 2925rga? - m^rfa\ - 22^rfa\ ^ 315r^^) 



J 



21r;^(af - r^)2(8af - 32r^a'j' 77r4af - 119r[;af -h 95r« 



(56) 



We also point out the consistency of our results with the numerical treatment of the Gregory- 
Lafiamme instability for rotating black string solutions [35] . There is a correspondence between the 
Gregory-Lafiamme instability of the black string and the existence of a thermodynamic negative 
mode of the black hole [36]. This correspondence implies that the negative eigenvalue is larger in 
magnitude if the threshold wavelength for the Gregory-Lafiamme instability is smaller, and this is 
indeed what is found when comparing our results with 
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C. Singly-spinning black ring 



The complexified singly-spinning black ring line element is [26J 



where 



F{x) 



dr - iCR^^dij 



+ 



R^F{x) 
{x - yf 



G{y) . dy2 dx2 ^G{x) . 



(57) 



F(0 = l + Ae, G(0 = (l-a(l + '^6, and C = ^ \{X - i^^j^^, 



(58) 



and the dimensionless parameters v and A lie in the range < z/ < A < 1. As it stands, the line 
element defined above has conical singularities aty = — 1, x = — 1 and x = 1. In order to remove 
the first two one must choose the periodicity of (p and ip to be 



A(/) = = 27r- 



1 



(59) 



The solution still has a conical singularity at x = 1, and is often referred to in the literature as the 
unbalanced black ring. To obtain the physically acceptable black ring one must remove the conical 
singularity at x = 1 by choosing A to be a function of u, leaving our solution dependent on two 
parameters {R,i'), 

2v 



A 



(60) 



1 + 1/2- 

The remaining parameters uniquely specify the mass and angular momentum of the black ring [32] . 
The singly-spinning Myers-Perry solution is obtained by taking the limit R ^ and A, z/ ^ 1 in 
the unbalanced solution, while keeping fixed a, tq given by 



2i?2 



and 



2i?^ 



X-u 



(61) 



" (l-i^)^ 
Following [26] , we note that the black ring instanton has another conical singularity located at the 
horizon y = —l/i', which is removed if we make the periodic identification (r, (p) = (t + /3, (p — ipQ), 
where p = 47ri?^Az^(l + X)/[VT^{1 + i^)] and n = ^/X^/[R^y X{1 + X)]. 

Again it is important to clearly state the range of coordinates in our manifold, which are given 
in the patch that we are interested in by r e [0,(3), (p^ip £ [0,27rVl - A/(l - u)), y £ 
and X £ [—1,1). We chose to determine I as a function of {R, u. A) and only latter imposing either 



relation (60 ) or the limit (61 ). The expressions for both 4>ab and X are complicated, but in the limit 



(61) T does reproduce Eq. (56|. Again we evaluated X using the relation (60), see Fig. (§, and we 
found that the black ring also exhibits an unstable behaviour against non-conformal perturbations. 
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Black Ring 
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FIG. 3: X is negative across the entire range. 



D. Six dimensional Myers-Perry 



In contrast to the previous cases, the moduh space of the six dimensional Myers-Perry solution 
is unbounded, in the sense that the black hole will exist for all values of one of the rotation 
parameters, if the other rotation parameter vanishes. 

It is a challenging problem to compute I for this black hole, since we are dealing with a co- 
dimension three manifold, so that finding a suitable line element for the computation is important. 
Here we will complexify the line element presented in [37], 

— ~ ~\~ z~i ~\~ 



+ 



R 
R 



(r^ _|_ y2^^j.,2 _|_ ^^2^ 

Y 



Y Z 

[dr - i{al - y^){ai - z^)d4>i - i{al - y^){al - z^)d4>2? 



[dr - i{ai + r^){al - z^)d</)i - i{4 + r^)(a^ - z^)d02 



2 I „2\/„ 2 



(r^ _|_ y2^^y2 _ ^2~) 

~ — wr^ — 2t[^'^ ~ ^("i + ~ - ^("2 + '^^)(«2 - y^)d(^2]^ 

where R, Y and Z are quartic polynomials given by 



(62) 



R(r) = {r^+al){r^+al)-4r, Y{y) = -{ai-y^){4-y^) and Z{z)=Y{z). (63) 
The coordinates (pi are related with the canonically defined coordinates {(pi,(p2) via 

and d>2 = n, .. . (64) 



af(af - a|) 
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I To - Six Dimensional Myers-Perry 
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FIG. 4: For the six dimensional Myers-Perry instanton, I is always negative. The dashed lines represent 
contours of constant Irp . 

The avoidance of conical singularities at the horizon, located at the largest root of R, implies 
the coordinate identifications (r, (pi, (1)2) = (r -|- (/>! — i/Jili, — iP^2), where /? = 47^rQr_^_/ R' {r_^_) 



and J7i = ai/{r^ + af). The region in moduli space in which the line element (62) represents a 
black hole, dependent on the zeros of R{r), is given by 

Wafalial - a\f - 4(af - ?>?>a\a\ - 33a|a? + a\)rl - TJrf < 0, (65) 

where the inequality is saturated for extremality and we have used the results compiled in [38]. 
This region of the moduli space obviously coincides with the region in which /3 is positive semi- 
definite. The range of the variables (r, r, y, z, ^2) is more involved: r G [0,/?), r G [r^,-|-oo), 
y £ [—b, b], z £ [b, a] and (pi £ [0, 27r), where b = min(|ai|, |a2|) and a = max(|ai|, |a2|). 

The numerical evaluation of I leads to a result very similar to the five dimensional case, see 
Fig. Q. T is negative for all values of (ai, 02) that have a positive semi-definite /?, i.e. black holes. 
The full expression considerably simplifies if we set both rotation parameters to zero, 

66 
'0 

which corresponds to the six dimensional Tangherlini solution [M]. Here, all extremal solutions 

have a non-divergent T, except when one of the rotation parameters vanishes and the other becomes 

infinite. Near the singular points, X is well approximated by 

7470^0 
280r, 



169rr 



^^-1^- (67) 



'0 
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E. Kerr-AdS black hole 



The Kerr-AdS instanton is given by 



^ (dr-4sin2M0)' + S2 



with 



dr2 d02 



S2 



adr 



— d(p 



(68) 



A = (r^ + a2)(l + 



ror, = + cos^ 6, Ag = 1 



cos 



and H = 1 



(69) 

where £ is the curvature radius of AdS and is related to the cosmological constant as £^ = 3/|A|. 
The Kerr solution is recovered in the usual limit i +oo. The line element (68) only represents 
a physical solution for \a\ < i, being singular in the limit |a| ^ ^ for which the 3-dimensional 
Einstein universe at infinity rotates at the speed of light 



The avoidance of a conical singularity at the horizon, located at r = r^, the largest root of A, 
requires, as in the Kerr case, the coordinate identification (r, (j)) = {t + (3 , (j) — z/3r2rot)- Here, 



47r(r^ + a? 



r+(l+a2^-2 + 3r2£-2 



a2r_|_2) 



and 



rot 



a(l -a2r2) 

J.2 _|_ g2 



(70) 



(71) 



The angular velocity Jlrot here is measured relative to a frame rotating at infinity. It is convenient 
to choose a coordinate system that is not rotating at infinity ((/;—>(/) — ia£~'^T), for which the 
angular coordinate identification depends instead on the angular velocity 



n 



(72) 



r2 + a2 

The reason for this is that Oj-ot is not the appropriate thermodynamic variable. It cannot be used 



to formulate the first law of thermodynamics [15j. The issue is irrelevant in this Section, as we 
are interested only in the coordinate independent quantity 2, but will be very important for the 
thermodynamic interpretation of the results. 

The expression for the TT perturbation (j)ab is remarkably simple, and is given by 



2S4 



2a^ sin2 9 2ia(r^ + o^) sin2 6 
1 + ^^T^ ^ ^ 






2iaE 



1 
-1 




S2H 






1 + 



2a^ sin2 6 
1]2 



(73) 
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FIG. 5: Phase diagram of the Kerr-AdS black hole. 



This perturbation is clearly traceless, and in the limit 



reduces to Eq. ( 49 ) . In the absence 



of rotation we know, through [11, that the negative mode ceases to exist exactly when the specific 
heat becomes positive. This transition occurs at = (3\/3/4)ro, the specific heat being positive 
for i < ic- This is similar to what happens when one considers Schwarzschild in a finite cavity 
|iO] . We expect that, in the absence of rotation, I will change sign for > £c, because (pat is not 
the negative eigenmode. In fact, Z changes sign at ij ~ 1.655ro, a value of i bigger, but close to 
ic- Classically, we also know that AdS black holes are unstable to superradiant modes whenever 

o£ > 1 m]. 

In Fig. ([5]), we represent a phase diagram for the Kerr-AdS black hole. In region I, I is positive, 
so we conclude nothing apart from the consistency with known results in the non-rotating limit 
|llj . The region of most interest is region II. Here the black holes are not superradiant, and may 
be classically stable, but I is still negative, meaning that black holes in this region are thermo- 
dynamically unstable. Region III corresponds to superradiant and thermodynamically unstable 
black holes. Finally, the white region corresponds to unphysical instantons, beyond extremality or 
the bound \a\ < L It is a curious fact that the line that divides non-superradiant black holes from 
superradiant black holes meets the contour line corresponding to X = at |a| = 
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V. THERMODYNAMIC INTERPRETATION 
A. Ensembles and stability 

The partition functions studied here have boundary conditions that impose periodicities both 
in time and in the rotation angles. This corresponds to fixing the temperature T = and the 
angular velocities 0,, i.e. to the grand-canonical ensemble. The canonical ensemble would require 
fixing the angular momenta J*, which includes specifying normal derivatives of the metric on the 
boundary. We expect a correspondence between the stability of the ensemble and the problem of 
non-conformal perturbations about the gravitational instanton. If there is a negative mode in the 
path integral, the partition function is ill-defined and the ensemble should not be stable. 

In this section, we recall the conditions for the thermodynamic stability of a system described by 
the grand-canonical ensemble. The laws of thermodynamics imply that, for a system in equilibrium 
at given temperature and angular velocities, any deviation from equilibrium, not necessarily small, 
obeys |l2] 

5E - T5S - VLi5f > 0. (74) 

Here, the angular velocities Jlj play the usual role of chemical potentials and the angular momenta 
J* are the respective charges (conserved in the canonical ensemble). Expanding 6E{S, ,P) to second 
order and using the first law of thermodynamics, 

dE = TdS + nidJ\ (75) 

we arrive at the conclusion that the stability relies only on the Hessian of the energy 



being positive definite. This Hessian matrix defines the so-called Weinhold metric ^3j. Its inverse 
is given by 

where G is the Gibbs free energy, 

G = E-TS-VLiJ\ (78) 

The positivity of the inverse Weinhold metric gives the same stability condition. Notice that the 
first law is equivalent to 

dG = -SdT - fdVLi, (79) 
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so that the coordinate systems and are related by = dE/dx^ and x^ = —dG/dy^. A 
third alternative is the positivity of the Ruppeiner metric 



g(R) = u^" = {E, r), (80) 

which is easily shown to be conformal to the Weinhold metric, 

d4 = - ^'^^f dn^du^ = /3 ^^^f dx^'dx" = I3dslr. (81) 

Let us now look at how the stability conditions relate to the usual linear response functions, 
like the specific heat. First, notice that 



dyij.dy^ 

where 



(82) 



is the specific heat at constant angular velocities (all ilj fixed). The isothermal differential moment 
of inertia tensor is 



There is also the vector 

where the second equality, given by the symmetry of the Hessian matrix, corresponds to a Maxwell 
relation. The rf dependence can be dealt with if we use the relation between the specific heat at 
constant angular velocities and the specific heat at constant angular momenta Cj, defined as 

C.-t[^), (86) 



which gives 



Hence 



Cn = Cj + T{e-%r,W. (87) 



q2q 

ds^ = — —dy^^dyu = pCjdyl + e'^iOiUjj, (88) 



where oJi = dyi + (e ^)ijr]Myo. 



From Eq. (88), we conclude that the thermodynamic stability is given by Cj and the spectrum 



of e*-' , the condition being 

Cj > and spec(e*^) > 0. (89) 
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B. Stability and negative modes 



We expect that the partition function correctly reproduces the thermodynamic features of a 
system. In particular, we expect that the partition function, which we associate with a system in 
equilibrium, has some sort of pathology if thermodynamic stability does not hold. In this section, 
we will extend a result by Reall |36] (based on previous contributions I45p on the relation 
between thermodynamic stability and the existence of a negative mode. He showed that a negative 
specific heat implied the existence of a negative mode in the action, in the case of the canonical 
ensemble. We wish to extend this argument to the grand-canonical ensemble. The converse result 
- to show that negative modes can only originate from thermodynamic instability - has not been 



proven, and we have found a counter-example (Section VC5|. 

The proof has two steps. The first step is to consider a path of geometries, intersecting a 
black hole saddle-point solution, on which the action functional has a certain dependence. That 
dependence implies an explicit relation between negative modes and thermodynamic stability. The 
second step is to show that the construction is possible, i.e. that such a path indeed exists. 

Consider paths of geometries, parametrised by variables for given temperature T and angular 



velocities Qj, for which the action takes the form 



/(z^r,%) 



(90) 



These paths intersect black hole solutions with temperature T and angular velocities which are 
given by = z!^{T,Qi). For the black hole geometries, / = PG, so that 

E{z^) = E{T,ni), s{z;:) = s{T,n,), j^«) = /(r,%). (91) 

The index in z^ has the same range as the index in = (T, i7j), the coordinates used in the 
last section to define the inverse Weinhold metric. The invertibility of zl^(y^) is assumed in the 
following. The argument works because the functions that only depend on can be thought of 
as depending on y^(z^). 

At saddle points of /, i.e. for black hole solutions. 



dl 
'dzi' 



dS 



^z^' V dy. 



dyu 



0, 



(92) 



which implies the first law of thermodynamics (75). The Hessian of / at these points takes the 
form 



dzt^dz' 



dT d^i 

dzt^ dzt^ 







- dT - 
















- 'dz^ - 



(93) 
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This result means that the negative modes of / given by this class of geometries have their origin 
in the non-positivity of the Weinhold metric, 



61 



5z^'5z'' = ^5sw^- (94) 



2 ^z^'^z 



2 

Notice that the path integral includes geometries that might not fall into this class. That is the 
reason why the proof of the converse result does not follow from the construction above. 



It would now be necessary to show that a path of geometries satisfying (92) is possible. Fortu- 
nately, this problem was addressed already by Brown, Martinez and York ^3] for the Kerr black 
hole, and the extension to higher dimensional cases should present no difficulty. They consider four 
dimensional axisymmetric geometries which satisfy the Einstein constraints and have the appro- 
priate boundary conditions on the wall of a finite cavity. In a such a cavity, it is possible to have 
thermodynamically stable black holes which are asymptotically flat (in this work, we consider the 
technically simpler boundary conditions imposed at infinity). The final result of the construction. 



their expression (43) for the gravitational action, exactly reproduces our assumption (92). 



Particular cases 



We have seen in Section V A that the linear response functions that characterise the stability are 
the specific heat at constant angular momenta and the moment of inertia tensor. In this section, 
we will calculate such quantities and interpret the results obtained with X. 

1. Kerr black hole 

Here the specific heat is given by 

3a* -I- 6r^a^ — 

where Cj is negative for small values of \a\, but becomes positive for \a\ > \/ 2\/3 — 3ro/2. The 
moment of inertia can be expressed as 

Cje = -^ (96) 

which is negative definite for black hole solutions. We conclude that, when Cj is negative, e is 
positive, and vice versa. From the grand-canonical point of view, the Kerr black hole is thermo- 
dynamically unstable, justifying why I was negative for all values of \a\ < ro/2. We have plotted 
both quantities in Fig. ([6|. Furthermore, in the absence of rotation, e is positive, meaning that the 
Schwarzschild black hole is stable against perturbations in the angular velocity. 
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FIG. 6: Graphic representation of Cj (solid line with axis on the left) and e (dashed line with axis on the 
right). 

2. Five dimensional Myers-Perry 



The specific heat at constant angular momenta and the moment of inertia tensor associated 
with the five dimensional Myers-Perry instanton are respectively given by 



Cj 



27r(x+ + af) + a^) {x^ - a{ 



alal[balal + 9(af + a|)x_^ + 20x2 Ix^^ + 3(af + al)x'^ 



(97) 



and 



(x, +a\){x, +a^)[3a|af +31 , (a|- 



aia2{x ,+a'^)(x ,+a1){a\a%-'ix'\_) 



2x 



-2x3^ 



2 — T~ 

3^ Gi-^ ^2 ~ 



aia2(x^+af )(x^+a|)(af a|-3a;^) 



{x^ +aj)(z^ +a|) [xj^ —Sa^x'^ +3afa2 (a|+2:^ 
2x,afa2—2x^ 



(98) 



We have plotted the region in parameter space for which e is not positive definite, and the region 
for which the specific heat at constant angular momenta is negative, see Fig. ([T]). As in the Kerr 
case, the two regions are complementary, meaning that the specific heat becomes positive when 
the moment of inertia ceases to be positive definite. This indicates that the Myers-Perry black 
hole is thermodynamically unstable in the grand-canonical ensemble, in agreement with the result 
obtained for I. 
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FIG. 7: The darker region indicates the region in the moduh space for which Cj is negative, and the hghter 
region indicates where e is not positive definite. 

3. Singly- spinning black ring 



The specific heat of the singly-spinning black ring is given by 



Cj 



(99) 



(l-Z^)(l + Z/2)3/2(2 + j,2)' 

where Cj is negative for small values of but becomes positive for v > 1/2. However, the moment 
of inertia, given by 

' = - (i-.)^aV.Y ■ 

is negative for all values of u. We have plotted both linear response functions in Fig. The 
fact that e is negative definite renders the black ring unstable in the grand-canonical ensemble, 
confirming the results obtained with Z. There is however a difference relative to all other black 
holes that we have analysed. In the ring case, the regions where Cj < and e < are not 
complementary. As a result, the specific heat at constant angular velocity, 

V2TT^R^i^^{iy{i^ + 2)-2) 



Cn 



(101) 



_ 1) + 1)3/2 ' 

does not have a definite sign, becoming positive for u > \/3 — 1. This means that Cq alone does 
not always identify the instability of the grand-canonical ensemble. We also point out that, in 
the extremal limit u = 1, both the specific heat and the moment of inertia are divergent. This 
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FIG. 8: Graphic representation of Cj (solid line with axis on the left) and e (dashed line with axis on the 
right). 

divergence is most likely related with the divergence of the energy and angular momentum at 
extremality. 



4-. Six dimensional Myers-Perry 



The specific heat at constant angular momenta is given by 



Cj 



27r(r2 + af)2(r2 + al)^[ajal - (a? + aj 



3r^ 



5af4 + Uajaliaj + 4)rl - {3af - Uajai + 3a|)r4 + 6(af + al)rl - 3r 



whereas the moment of inertia tensor is e = (r^ + )(r^ + a^e, where 



3rl - (a2+6a2)r4 +3iaf+2a^al )r2 +3afal 



2r. [3r4 -(a2+a|)r2 



11012 [^r-^ -(af +£12)''^ -afail 



aia2[7r4 -(a2+a|)r2 
' r, [3r4-(af+ai)r2 



_[3rl-{al+al)rl-a 



if ail 



2r-+[3r4-(af+ai)r2-afai] 



(102) 



(103) 



We have again plotted Cj as a function of the normalised rotation parameters ai/ro, see Fig. (|9]). 
The similarities with the five dimensional Myers- Perry solution are evident. The region in which 
Cj is negative is complementary to the region where one of the eigenvalues of e is negative. We are 
again led to conclude that the six dimensional Myers-Perry solution is thermodynamically unstable 
in the grand-canonical ensemble, in agreement with the results obtained for X. 
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six Dimensional Myers-Perry 




FIG. 9: The darker region indicates the region in the moduh space for which Cj is negative. The region in 
which e is not positive definite includes the hghter region. 

5. Kerr-AdS black hole 



The Kerr-AdS case is richer but also more subtle. We can calculate the specific heat and moment 
of inertia using the expressions in [15j, for which we obtain 

27r^4(Q2 ^ rl)\rl - f)a^ + r\{e + Sr^ )] 



Cj 



(a2 -£2)[(3a4 + 6r2a2 



+ (a6 + 13r2 a4 + 2'ir%a^ + Sr^ )^2 + a^r\{a^ + 3r2_)2] 



and 



7r£8(a2 + r2 )3[(r2 _ ^2)^2 + ^2 (3^2 + ^2)] 



(104) 



(105) 



r+(a2 - l^Y\{r\ + ^2)a2 _ r2_(3r2 - £2)] • 

The parameter region for which the last expression is negative is composed, as in the previous 
Kerr and Myers-Perry cases, of two complementary regions. The specific heat at constant angular 



momentum is negative in one of them - regions I and IV in Fig. (10) - and the moment of inertia 



is negative in the other - regions II and V in Fig. (10). However, there is a region where both are 



positive - regions III and VI in Fig. (10) - which contains, in the non rotation limit a = 0, the 



known stable phase of large Schwarzschild-AdS black holes. 



The areas IV, V and VI in Fig. (10) identify the gap between a negative mode, found in 



Section IV E and the line of critical stability, found using the thermodynamics formulae. On the 
left, we have a strip in which we expect thermodynamic instability, based on the linear response 
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FIG. 10: Phase diagram of the Kerr-AdS black hole. In region III, the linear response functions Cj and e 
are positive and we find I > 0. In regions II and V, e < and C,j > 0, while, in regions I and IV, Cj < 
and e > 0. We find a negative mode, i.e. X < 0, in regions I, II and VI. The thin region VI, amplified in the 
top-right corner, represents a puzzling feature since we find a negative mode when thermodynamic stability 
was expected. 

functions Cj (in region IV) and e (in region V), but find no negative mode. This is not surprising 
because our perturbation is not a negative eigenmode and we do not expect it to identify the 
entire region of instabihty, which in the non-rotation Umit a = persists until the specific heat 
becomes positive. There is however a very thin strip, region VI, amplified on the top-right corner 



of Fig. (10), for which we detect a negative mode when no instability was expected in terms of 



the thermodynamic response functions. We find X < but the Weinhold metric indicates stability. 
This is a puzzling feature as it seems to contradict our expectations that the partition function 
exactly reproduces the known thermodynamics. 

The statement of the puzzle above immediately indicates a possible explanation for this small 
discrepancy. It is well known that the partition function, canonical or grand-canonical, which 
accounts for statistical physics, does not always reproduce thermodynamics. The 'thermodynamic 
limit' arises only as a macroscopic approximate description. In the grand-canonical ensemble, the 
partition function reads 



Z{T,n) = j dEdJuj{E,J) e 



-l3{E-nj) 



(106) 



where uj{E, J) is the density of states with given energy E and angular momentum J. The variance 
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of the statistical distribution is given by 

{E^)-{E)' _ fd{E)\ 

jJ') - jJ)' _ T (d{J)\ 

These variances are not small in most of the parameter space and diverge along the critical line 
of stability if we take {E) and (J) to be the thermodynamic quantities E and J used before. The 
derivative factor on the RHS of ( |108 ) is then the moment of inertia e, and the quantity held fixed 



in the derivative on the RHS of (107) is related to the so-called fugacity in the statistical mechanics 



of systems of particles. We saw that the RHS of (108) changes sign across the divergence on the 



line of critical stability, but the RHS of ( 107 ) is always positive in spite of the divergence. 

In the region close to a phase transition, a careful study of the contribution of fluctuations is 
essential. If the thermodynamic description breaks down near a critical point, the precise deter- 
mination of quantities like the critical temperature Tc requires the contribution of higher order 
corrections to the partition function. A well understood example is the two-dimensional zero-field 
Ising model, which can be solved exactly. The exact critical temperature corresponds to almost half 
of the naive leading order result, and the first correction is already very important (see Fig. (9.5) 
of [l6]; the van der Waals gas is another good example described in this reference). In our case, the 
one loop contribution is the first quantum correction and thus it is not surprising that the critical 
line, as predicted by the zero-th order quantities, will suffer a small correction. 

We must point out that a linear response function, namely the specific heat at constant charge, 
also diverges in other known cases and the negative mode still disappears when expected. This hap- 
pens for the Schwarzschild-AdS black hole, for which the specific heat diverges at £c = (3\/3/4)?"o, 
where the negative mode has numerically been shown to disappear . It also occurs in the case of 
the Reissner-Nordstrom black hole, for which the specific heat at constant electromagnetic charge 
Cq diverges at Q = where the negative mode has analytically been shown to disappear 

|47j (the latter work focuses on the canonical ensemble, so that both the stability and the validity 
of the thermodynamic limit are given by Cq). In the Kerr-AdS case studied here, we are looking 
at the grand-canonical ensemble. The specific heat at constant angular momentum Cj is finite at 



the transition, which is signalled instead by the quantities on the RHS of (107) and (108). Since 
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or, taking (E) and (J) to be the thermodynamic quantities E and J used before, 

-c,^,.m\ ,no) 



Now, Cj is divergent and e = along the hne separating the regions I-IV and II-V in Fig. ( 10 ) . This 
corresponds to the transition in the canonical ensemble. In the grand-canonical ensemble, this line 
has no special interest and the LHS of the latter expression is still finite, because the divergence in 
{dE / dJ)T cancels the divergence in Cj. On the other hand, the line of critical stability, separating 
regions III- VI from the rest, is marked by a divergence in e that also causes the divergence on the 
LHS, while the other quantities are finite. The exception is the single point for which a = 0, the 
Schwarzschild-AdS case. There, the transition is signalled by the divergence of the specific heat 
since the moment of inertia is actually finite for a = 0, e = rJ|„/2. 

We conclude that the phase transition in the grand-canonical ensemble can be traced back to 
the divergence of the moment of inertia e only. We speculate that this transition may have different 
quantum properties such that the critical stability line is corrected by the one loop contribution, 
while the same does not happen in the canonical case with the specific heat at constant charge, in 
the cases known. 

To complete the discussion of the transitions, let us consider the critical exponents. Notice 
that these too can be corrected by higher order contributions to the partition function, but we 
are only able to study them here for the zero-th order semiclassical result corresponding to our 
thermodynamics formulae. 

We analysed, in the canonical ensemble (charges J or Q fixed), the cases of Kerr black hole, 
with or without cosmological constant, and of the Reissner-Nordstrom black hole. We use the 
definitions C ~ ib|T — Tc|~", for the phase T > Tc, and C ~ ib|T — Tc|~"', for the phase T < Tc- 
We will actually find that the critical temperature is a local maximum or minimum as a function 
of the energy E^ for fixed charges J or Q. This means that only one of the exponents a and a' is 
relevant. In a physical situation, on one side of the phase transition there is the stable black hole, 
while on the other side there is a different stable configuration (eg. spacetime with radiation j40j ) 
the specific heat of which would contribute with the other critical exponent. In the asymptotically 
flat cases, Tc is a maximum of the temperature as a function of the energy E, for fixed J or Q, 
which means that the relevant exponent is a' . In the Schwarzschild-AdS case, Tc is a minimum of 
the temperature as a function of the energy E, so that the relevant exponent is a. In the Kerr-AdS 



case, both situations can occur depending on the parameter region, as can be seen in Fig. (11). 
Since, in all cases, the temperature varies quadratically with the energy around the critical point. 
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FIG. 11: In the Kerr-AdS case, up to three black holes can exist with a given temperature T, for a certain 
range of a fixed angular momentum J. In this example, J — 1 and t — 10. In fact, a fourth solution exists 
with small values of the energy E, but is beyond the limit \a\ — I. 

T — Tc ^ ±{E — Ec)^^'^, the critical exponents have the value 1/2. 

In the grand-canonical ensemble, the charge is not fixed. The interesting cases are the asymp- 
totically AdS black holes since there is no thermodynamic stability for the asymptotically flat black 
holes in this ensemble. In the Schwarzschild-AdS case 17 = 0, the specific heat is the same as in the 
canonical case, diverging as described above. We have also mentioned that the moment of inertia, 
given by e = is finite. When $7 7^ 0, there is a sharp difference as e becomes divergent over 



the line of critical stability. Let us recall equation (87) in four dimensions, 

Cn = Cj + -^. 



fill) 



The quantity r] diverges as e and Cq, while Cj is finite, along the critical line. Again the divergences 
have a critical exponent 1/2, at fixed 0, (and the same dependence on 0, — 0,c for fixed T). 



VI. CONCLUSIONS 



In this paper, we studied the thermodynamic instabilities of several rotating black holes. We 
analysed their grand-canonical partition functions and found that they possess negative modes 
at the perturbative level. An argument connecting this pathology to the lack of thermodynamic 
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stability is extended here to the grand-canonical ensemble. 

The method we applied to look for negative modes avoids the complications due to the lack of 
symmetry of these solutions. Instead of addressing the full problem, we construct a non-conformal 
probe perturbation which lies along the stationary phase path. If the probe perturbation decreases 
the Euclidean action, then a negative mode exists. We also clarified the issue of using Euclidean 
methods for rotating spacetimes, leading to complex (quasi-Euclidean) instantons which should 
present no difficulty of principle. 

We found our results for the negative modes to be consistent with the standard conditions of 
thermodynamical stability in the grand-canonical ensemble, with the exception of the Kerr-AdS 
black hole. In the latter case, a small parameter region near the line of critical stability presents 
a negative mode when the zero-th order conditions of stability are satisfied. This feature suggests 
that the thermodynamic limit is ill-defined near the critical line in this case, so that the location 
of the phase transition is corrected by quantum perturbations. 

Possible future work linked to these results includes: (i) studying the partition function of a dual 
conformal field theory and comparing with the gravitational case; (ii) studying the full negative 
modes problem numerically, with special focus on the phase transition of the asymptotically AdS 
black holes; (iii) formulating the canonical ensemble perturbations problem, which would require 
fixing the angular momenta J*, i.e. writing the second-order action in terms of perturbations that 
would specify derivatives of the metric normal to the boundary. 
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